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We present numerical Monte Carlo results for the stationary state properties of KPZ type growth 
in two dimensional surfaces, by evaluating the finite size scaling (FSS) behaviour of the 2nd and 4th 
moments, W2 and W4, and the skewness, W3, in the Kim-Kosterlitz (KK) and BCSOS model. Our 
results agree with the stationary state proposed by Lassig. The roughness exponents W n ~ L an 
obey power counting, a n = na, and the amplitude ratio's of the moments are universal. They have 
the same values in both models: W 3 /Wa B = -0.27 (f) and W 4 /W% = +3.15 (2). Unlike in one 
dimension, the stationary state skewness is not tunable, but a universal property of the stationary 
state distribution. The FSS corrections to scaling in the KK model are weak and a converges well 
to the Kim-Kosterlitz-Lassig value a = |. The FSS corrections to scaling in the BCSOS model are 
strong. Naive extrapolations yield an smaller value, a ~ 0.38(1), but are still consistent with a = | 
if the leading irrelevant corrections to FSS scaling exponent is of order y ir ~ —0.6 (2). 

PACS numbers: 02.50.Ey, 05.40+j, 68.35.Fx 



I. INTRODUCTION 

KPZ type growth is one of the generic dynamic pro- 
cesses describing growth of crystal surfaces. It is named 
after the Langevin equation introduced by Karder, Parisi, 
and Zhang about a decade ago 
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with 77 uncorrelated noise 
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Numerous microscopic models on the master equation 
level have been studied numerically as well and are con- 
firmed to be in the KPZ universality class || . However, 
many properties of this process are still in question, in- 
cluding basic aspects, like the precise values of the scaling 
critical exponents, and the detailed structure of the sta- 
tionary growing state. Part of the problem is the absence 
of an obvious mean field theory. The linear, integrable 
diffusion, part of the KPZ equation (A = 0) does not play 
the role of mean field theory fixed point in high enough 
dimensions (D). The KPZ behaviour is governed by a 
strong coupling fixed point for all D, and thus evades 
perturbative renormalization treatments . 

It is widely accepted that the dynamic exponent z 
and the stationary state roughness exponent a obey the 
equality a + z = 2 in all D These critical exponents 

specify how time and height rescale under a renormaliza- 
tion transformation: x — > bx, t — > b z t, and h — > b a h. 
The exponent identity states that under renormalization 
the amplitude of the non-linear term in the KPZ equa- 
tion does not change. A is a so-called redundant scaling 
field. It plays a role similar to lattice anisotropy in equi- 
librium phase transitions. Increasing A simply speeds-up 
the process, and scaling amplitudes are proportional to A. 



This exponent equality links the dynamic scaling to the 
stationary state scaling. Therefore the focus has shifted 
recently to the structure of the stationary state. 

The stationary growing state is trivial in one dimen- 
sion {ID). It is the Gaussian distribution. The up and 
down steps along the surface are uncorrelated beyond a 
definite correlation length. This implies a — j (the ran- 
dom walk value), and from the above exponent identity 
it follows that z=\. This behaviour is well established, 
not only by numerical studies |^-||] , but also analytically. 
The ID body centered solid on solid (BCSOS) growth 
model is exactly soluble Jt],||]. It's master equation is a 
special case of the 2D equilibrium 6- vertex model. In the 
latter representation, KPZ scaling describes facet ridge 
end-points of equilibrium crystal shapes ||. ID KPZ 
growth is equivalent also to asymmetric exclusion hop- 
ping processes [[To). Moreover, the exact stationary state 
of the Langevin equation itself is known in ID and is 
indeed the Gaussian distribution [^J). 

The stationary state is not simple in D > 1. The sta- 
tionary state roughness exponent a takes a non-trivial 
value and is actually not very well known numerically. 
For example, in 2D the reported values vary between 
a = 0.37 — 0.4 P,^,[TT|-^3t . Lassig made an important 
analytical break through last year |L4[ . He proposed the 
likely structure of the stationary state by studying the 
operator product expansion of the (height variable) cor- 
relation functions. Under the assumption that the alge- 
bra closes and contains only one scaling field operator, 
Lassig obtained a quantization condition for the expo- 
nents. One of these solutions, a = |, is close to the 
above 2D numerical values. The moments 



W n = ((hi-h) n ) 



(3) 



of Lassig's stationary state distribution obey power 
counting, i.e., the exponents in the scaling relations 



W n (N- L ,t- L ) =b a "W n (bN- L ,b z t- L ) 



(4) 



are related as a n = na. The distribution lacks multi 
scaling. Moreover, for the closure of the algebra it is 
important that the stationary state is skewed. The odd 
moments, in particular the third one must be non-zero. 

In this paper we report a detailed numerical study 
of the stationary state in 2D for the Kim-Kosterlitz 
(KK) fflj and the BCSOS model. We determine the 
finite size scaling behaviour of the second, third, and 
fourth moments. 

In section 2 we review the properties of stationary state 
skewness in ID KPZ growth, and list the possibilities 
in higher dimensions. Section 3 contains our numeri- 
cal results for the KK and BCSOS model. The third 
moment is indeed non-zero and the second, third, and 
fourth moments obey indeed power counting. We find 
strong evidence that the amplitude ratio's of the mo- 
ments are universal, R 3 = W^/W^ 5 = -0.27(1) and 
i? 4 = W±/W% = +3.15(2). The amount of stationary 
state skewness is the same in these two models. In ID 
skewness is tunable, but appears not to be so in 2D. The 
values of these amplitude ratio's are the major new re- 
sults of this paper. 

In section 3 we address also whether the difference in 
the previously reported values of a for the KK and BC- 
SOS models can be attributed to finite size scaling (FSS) 
type corrections. The FSS corrections to scaling expo- 
nent is large, but in the same range of values suggested 
by naive power counting. 

To check more directly whether skewness is tunable or 
not, we introduce a temperature type parameter K in 
the BCSOS model. In the KK model it is known that A 
changes sign with K jl2| . In section 4 we give an intuitive 
explanation for why this happens. (It is related to pre- 
roughening phenomena in equilibrium surfaces.) In the 
BCSOS model A does not change sign. In section 5 we 
present MC data for the K dependence of the roughness 
exponents and the amplitude ratio's. Both show some 
systematic drift, but much smaller than expected if they 
would vary with K. 



II. STATIONARY STATE SKEWNESS 

This study was actually not motivated by Lassig's re- 
cent results. It was conceived as a generalization of an 
earlier study of stationary state skewness in ID The 
KK model is a special point in the restricted solid-on- 
solid (RSOS) model. We varied its adsorption and evap- 
oration probabilities in the ID model by making them 
dependent on the local nearest neighbour heights. That 
lead to 5 independent parameters. Surprisingly we were 
able to construct the exact stationary state in a 4 dimen- 
sional subspace. Its structure is simple. The steps in the 
interface are completely uncorrelated. Only the step den- 
sity varies. This state has zero skewness. It is a Gaussian 
and has particle-hole symmetry. Outside this exact solu- 
ble subspace the stationary state is skewed. This means 



that in general KPZ type growth in ID has a non-zero 
third moment W3 in its stationary state. For example, 
the KK point lies outside the non-skewed subspace. On 
the other hand, the stationary states in the exactly sol- 
uble BCSOS model and also the Langevin equation it- 
self are non-skewed. Stationary state skewness is distinct 
from temporal skewness. Most initial states develop into 
skewed structures at intermediate time scales (temporal 
skewness) even if the stationary state is not skewed [[l6| . 

So in ID, KPZ type stationary states are typically 
skewed, and its amplitude is tunable. This raises the 
immediate question whether skewness affects the scaling 
exponents. Let's refer to the operators leading out of the 
non-skewed subspace as O sk and to their conjugate cou- 
pling constants as u sk . In ID, a KPZ type fixed point 
with zero-skewness exists. The question is whether this 
fixed point is stable; whether O sk is a relevant, marginal, 
redundant, or an irrelevant operator. Suppose the non- 
skewed KPZ fixed point is stable. The moments should 
scale then with system size as 

W n (N-\u sk ) = b a "W n {bN~\b^u sk ) (5) 

with q„ = na and y sk < 0. All even moments scale as 

W n ~ AN a " + .... (6) 

with universal amplitude ratio's, R n — Wn/W^ 2 ■ All 
odd moments scale as 

W n (N~\u sk ) = N a -T n (N y ° k u sk ) (7) 
= N a - [T n (0) + ^(0)Nv«*u sk + ...] 

~ Usk N a -+y^ 

with J- n (0) = because the fixed point has no skewness. 
The odd amplitude ratio's are proportional to u sk , i.e., 
the skewness varies continuously. 

Numerical (transfer matrix finite size scaling) results 
confirmed that the non-skewed KPZ fixed point is stable 
in ID. We found y sk ~ — 1. Moreover, the amplitude of 
W3 (at e.g., the KK point) is indeed roughly proportional 
to the skewness coupling constant u sk , in accordance with 
eq-0 

Skewness is negative at the ID KK point. On average 
hill tops are wider (flatter, less sharp) than valley bot- 
toms. Such a statement is meaningless without a spec- 
ification of a cut-off. The definition of what constitutes 
a mountain and what represents a local hump depends 
on the length scale at which the surface configuration 
is being viewed. (Humans do not interpret every grain 
of sand as a hill.) Skewness is a scale dependent prop- 
erty. The asymptotic scaling of the moments tells us how 
asymmetric the hills and valleys are in the large length 
scale limit. In ID this skewness persists all the way to 
microscopic length scales. We calculated the surplus of 
sharp valley bottoms over sharp hill tops at the micro- 
scopic (the grains of sand) level. This expectation value 
has the same sign as the macroscopic skewness and is 
also roughly proportional to u sk . This invariance of the 
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surface structure over all length scales is related to the 
rather trivial nature of the fixed point stationary state 
(Gaussian distribution). 

The tunability of the skewness at the microscopic 
length scale is easy to understand. Consider the ID 
RSOS model with deposition only, with three parame- 
ters: ph, p Sl and p v Jlq| . The density of local sharp hill 
tops is set by the deposition probability p^ of particles on 
local flat surface segments. The sharpness of these local 
hilltops is set by the rate at which they broaden, i.e., the 
probability with which particles adhere to existing steps, 
p s . The density of local sharp valleys is set by the rate 
at which single particle puddles fill-up, p v compared to 
the rate at which they are created, 2p s . These processes 
balance exactly inside the subspace where the station- 
ary state is trivial and has zero skewness. At the KK 
point, ph = p s = Pv = l-i the balance is imperfect and 
the dynamics creates a backlog of "to-be-filled-up" local 
valleys. Newly created local hill tops broaden readily, 
and therefore are flatter [^5) . 

Is the skewness tunable in D > 1 as well? Is there 
maybe a line of KPZ type fixed points with continuously 
varying skewness? A varying exponent a would explain 
the current numerical spread in its value. Is it possible 
to change the sign of skewness without changing the sign 
of A in eq. ([[])? Or. is there one single KPZ stationary 
state fixed point distribution, with a universal amount 
(and sign) of skewness? In that case we need to explain 
the numerical spread in a' s in terms of strong FSS cor- 
rections. These are the issues we address in this paper. 

III. SCALING OF THE STATIONARY STATE 
MOMENTS 

We perform a systematic numerical study of the sta- 
tionary state properties in the 2D KK-model and the 
BCSOS model. In both cases we allow only particle de- 
position (no desorption). Consider a 2D square lattice 
with an height variable h(r) = 0, ±1,±2, ... at each lat- 
tice site. We apply periodic boundary conditions. In 
the KK-model, nearest neighbour columns are allowed 
to differ by at most one unit, 5 h — 0, ±1. Choose a site 
at random, and deposit a particle, h — > h+ 1, with prob- 
ability p = 1 unless such a move would violate the above 
restriction. 

In the BCSOS model the square lattice is divided into 
two sublattices. The height variables are restricted to be 
even on one of them, h(r) = 0, ±2, and to be odd on 
the other h(r) — ±1, ±3, such that nearest neighbour 
columns always differ in height by one Sh = ±1. Choose 
at random a site, and deposit a particle with probabil- 
ity p = 1 if the move does not violate the Sh = ±1 
constraint. In section ^ we consider also the generaliza- 
tion where the BCSOS deposition probabilities vary with 
the local configuration by means of a temperature type 
parameter, p = min(I, exp(— AE)). AE is the energy 
change if the adsorption would take place. The energy 



E({hi}) = E \ K V* - h >f ( g ) 

has a tunable parameter K and the summation runs over 
all next nearest neighbours. This rule is identical to 
that in standard Metropolis MC simulations, except that 
desorption is forbidden. The latter breaks detailed bal- 
ance and leads to the non-equilibrium growing stationary 
state. 

We determine the second, third, and fourth moments 
of the stationary states. The MC averages in this section 
involve ~ 2 10 6 MC steps, after ~ 4 10 3 initial MC config- 
urations, to allow the surface to reach its stationary state. 
The square lattice size 1? varies between 12 < L < 128. 

First consider the KK model. Fig. |] shows the sec- 
ond moment W 2 . It scales as W 2 — AL° 2 . The slope 
of the log-log plot gives the exponent a 2 - It is a mis- 
take to apply a least-square type fit to the slope at 
large L. One should determine the slope at various 
system size intervals and perform a finite size scaling 
(FSS) analysis. Fig. || represents such an analysis. It 
shows FSS estimates for a 2 from the same data, de- 
fined as a 2 (L) = \og[W 2 {L 2 )/W2{L l )]/\og[L 2 /L l ] with 
L 2 ~ 1.2L\ and L = \{L X + L 2 ). Such a FSS scaling 
analysis is only barely feasible at large system sizes due 
to the intrinsic MC noise. Fig. ||is indeed rather noisy at 
large L. The MC scatter increases with system size, since 
the stationary state is intrinsically critical and therefore 
subject to critical slowing down. We opted for running 
many lattice sizes instead of fewer but longer MC runs. 
The solid line is obtained by averaging the a 2 (L) locally, 
over L — 7 < L < L + 7. The data in Fig. || converge 
by eye to a 2 = 0.80(2), consistent with the value a = | 
proposed by Kim and Kosterlitz from their earlier nu- 
merical results , 0| , and with Lassig's stationary 
state. The FSS corrections to scaling in the KK model 
are small compared to the MC noise. 

Figs. H and [| show the same FSS analysis for the third 
and fourth moment exponents 03 and 04. They con- 
verge by eye to 03 = 1.20(4) and 0:4 = 1.60(5). This 
agrees with Lassig's stationary state. In particular, the 
surface is skewed, and power counting, a n = not, is sat- 
isfied within the numerical accuracy. The corrections to 
FSS are again small compared to the MC noise. 

Figs. H and || show the amplitude ratio's (the circles) 
of the third and fourth moment compared to the sec- 
ond one, R n = W n /(W 2 ) n/2 . The fact that these ratio's 
convergence gives additional evidence of the validity of 
power counting. Actually, they do so much smoother 
than the exponents a n , suggesting that the MC fluctua- 
tions tend to preserve the power counting property better 
than the precise value of a. The amplitude ratio's con- 
verge smoothly to R 3 = -0.27(1) and R A = +3.15(2). 
Skewness is negative, R3 < 0. It is probably wishful 
thinking to guess that R4 — it. (R4 = 3 in Gaussian 
distributions.) 
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Consider the BCSOS model at K = 0. Figs. @, |, 
and [I show the same type of FSS estimates for the expo- 
nents a n . The finite size corrections to scaling are much 
larger than in the KK model. The data converge by eye 
systematically to smaller values than in the KK model: 
a 2 = 0.77(2), a 3 = 1.16(3), and a 4 = 1.54(4). This 
is consistent with the value for a 2 reported in the lit- 
erature p jl3|. These a n 's are mutually consistent with 
power counting. Power counting is again more stable 
than the values of the exponents a n . Most importantly, 
the amplitude ratio's in Figs. |5| and || (the diamonds) 
converge to the same values as in the KK model (the 
circles) . 

Is a really smaller than in the KK model and differ- 
ent from the Kim-Kosterlitz-Lassig (KKL) value? The 
finite size scaling corrections in the BCSOS model are 
several orders of magnitude bigger than in the KK model. 
Is it believable that these apparent differences are due 
to corrections to scaling only? The above FSS extrap- 
olation "by eye" presumes implicitly that all a n con- 
verge approximately linearly in L . This looks rea- 
sonable from the data, but is too restrictive. Correc- 
tions to scaling originate from so-called irrelevant scaling 
fields. The corrections to scaling exponents yi r < in 

W n ~ A n N an [1 + B n L Vir + ] are universal properties 

of the stationary state fixed point. The amplitudes B n 
are not universal. They depend on the "distance" of the 
model to the fixed point and the quantity we are looking 
at. Assume that one correction to scaling term domi- 
nates, i.e., that all other operators scale with much more 
negative values of yi r . The same exponent yi r should 
then appear in all moments. The only exception is that 
in specific quantities the leading term might have zero 
amplitude by symmetry. For example, all even moments 
might show a different leading exponent yi r than all odd 
ones. This actually happens here. 

Suppose we force our BCSOS data to converge to 
a = f. We made plots of W n /L na with a = § ver- 
sus L Vir for a range of values of yi r . This should be a 
straight line at the proper yi r . From this we estimate 
that yi r ~ —0.6 (2) for the second and fourth moments, 
and yi r ~ — 1.7 (3) for the third moment. These straight 
line fits are satisfactory stable. So the KKL value for a 
is within the realm of possibilities for the BCSOS model. 
Still, it remains a leap of faith, because the curves in 
Figs. |7|, [|, and ^|, are bend to the limit. 

It would be much more convincing if the corrections 
to scaling exponents where known analytically, and/or 
take simple values. At the core of Lassig's result is the 
assumption that the operator content of the system is 
simply. Therefore one would expect that the irrelevant 
operators have rather trivial critical dimensions, like in- 
tegers, multiples of a, and combinations of both. The 
above numerical values of yi r do not look that simple, but 
are of the same order of magnitude as we might expect. 
Simple minded power counting in Eq.(|l|) with a = | 
and z + a = 2, suggests that the corrections to scaling 
are strong; that the curvature operator d 2 h/dx 2 is irrel- 



evant, but not by much, y v = —a, and that the leading 
skewness operator (d 2 h/dx 2 } scales with y s k = —2. 



IV. TEMPERATURE DEPENDENT 
TRANSITION PROBABILITIES 

In section [v| we vary the temperature type parameter 
K in the BCSOS model, to study the universality of a 
and skewness issues raised in section III in more detail. 



We do this only for the BCSOS model, not the RSOS 
model. (The KK model is a special point in the latter.) 
From earlier studies it is known that in the RSOS model 
the KPZ non-linear term A changes sign with K [ Q . 
This creates strong Edwards- Wilkenson (EW) type cor- 
rections to scaling and will obscure the skewness prop- 
erty, because the EW stationary state at A = is non- 
skewed (the Gaussian distribution). Such a change in A 
does not take place in the BCSOS model. 

The following connection with preroughening phenom- 
ena in equilibrium crystal surface explains why A changes 
sign in the RSOS model and not in the BCSOS model. 
Imagine a two dimensional phase diagram, with K ~ T^ 1 
the temperature like parameter, and a parameter s rep- 
resenting the asymmetry between particle deposition and 
evaporation, s = corresponds to MC equilibrium type 
dynamics and s = 1 to the above pure deposition model 
without evaporation. The equilibrium surface undergoes 
a roughening transition. The rough phase is the EW sta- 
tionary state. The scaling properties of the equilibrium 
(stationary) state are described by the Gibbs distribution 
of the sine-Gordon model 



E = dx dy 



K 



— (V/i) 2 + u 2 cos(27i7i) + u 4 cos(47r/i) 



(9) 



It is known that u 2 varies with temperature and changes 
sign inside the equilibrium rough phase in the RSOS 
model just above the roughening transition. This follows 
from an exact duality transformation JT^]. The location 
of this u 2 = point at s = agrees qualitatively with 
the numerical value of the A = point at s = 1. Those 
values do not have to be identical. They only need to 
be of the same order of magnitude. The following argu- 
ments suggest that a line of A = points emerges from 
the u 2 = point at s = into the s direction. 

u 2 is negative at the high temperature side of the 
u 2 = line. There, the rough stationary state surface 
takes locally the so-called disordered flat type structure 
with alternating up and down steps and local half integer 
surface height |l^] . This is the same surface structure as 
in the BCSOS model (but not close-packed with dh = ±1 
kinks). The non-linear term in the KPZ equation con- 
trols the local growth velocity at sloped sections of the 
surface. Growth at slopes is suppressed in BCSOS type 
rough structures, and therefore A < 0. At the opposite, 



4 



low temperature ui > side of the U2 — point, the lo- 
cal rough RSOS surface is smooth, and has on the local 
level integer average surface heights. In such structures, 
growth at sloped parts of the surface is enhanced and 
therefore A > 0. 

The location of ui = can be controlled by introduc- 
ing further neighbour interactions. This point transforms 
into the preroughening transition point when it moves 
below the roughening temperature |lq| . (But only at 
s ^ 0, because the driven non-equilibrium surface is al- 
ways rough). 

Such a change in A does not take place in the BC- 
SOS model. Its equilibrium stationary state (at s = 0) 
is exactly known for all K from the exact solution of the 
U2 is negative at all values of K 
equilibrium surface. Therefore there is no 



6-vertex model |19 
in the s 



reason to expect a change of A as function of K in the 
pure deposition model at s — 1. This makes the BCSOS 
model a suitable testing ground for the tunability of 2D 
skewness and the universality of a. 



V. TEMPERATURE DEPENDENT DEPOSITION 
RATES IN THE BCSOS MODEL 



In this section we study the universality of the sta- 
tionary state roughness exponent a and the skewness by 
varying the deposition probabilities in the BCSOS model. 
This should clarify whether the skewness is truly univer- 
sal or a tunable parameter. We can control the micro- 
scopic particle-hole asymmetry explicitly. 

Fig. " 



10 shows the variation of the second moment ex- 



ponent Qt2{L) with the temperature parameter K. This 
plot is more qualitative than the ones in section [II . The 



system sizes are smaller (up to L=80) and the MC runs 
an order of magnitude shorter (up to 10 5 MC steps, after 
2500 initial MC steps). The curves have a definite slope 
at small system sizes, but these disappear with system 
size. At temperatures beyond K ~ 1 the surface be- 
comes very flat and inactive. Compared to K ~ the 
system size is effectively much smaller and the MC runs 
effectively much shorter. This explains the decay in the 
approximants at large K. At the opposite side, beyond 
K ~ — 1, the surface becomes quite facetted at short 
distances and the dynamics slows down again. Facetted 
structures have a = 1. This explains why the aa(i) 
curves drift upward on the left hand side in Fig. a 3 
and ot\ vary also only weakly with K. The amplitude 
ratios R 3 = W^/W^ 5 and i? 4 = W4/W2 do not vary 
significantly with K either. 

We performed quantitative MC runs at K = 
±0.25, ±0.1 for system sizes up to L = 60 (with 2 10 6 
MC step runs after 4 10 3 initial MC steps). Figs. ||, 
O, and [l3] show the a n (L) approximants as function of 
l/L at K = ±0.25. The K = data is included as ref- 
erence. The corrections to scaling in the even moments 
increase with K , and are an order of magnitude smaller 



at K — —0.25. The three curves are consistent with con- 
vergence towards the same values for 0.2 and a±. Naively 
these values point to an a smaller than | , but still consis- 
tent with the KKL value if the crossover scaling exponent 
is large (see section III). The corrections to scaling in a 3 
are less clear cut. At first glance the three curves seem 
more or less parallel (which would suggest a continuous 
variation in 03 with K) 7 but they actually converge to 
indistinguishable values for given the error bar from 
the numerical MC noise. 

3/2 

The amplitude ratio's R 3 — W3/W0/ and R4 — 
W4/WI are shown in Figs. [l4| and [l5]. Like before these 
amplitude ratio's are more stable than the numerical val- 
ues of a n . The data confirm universal ^-independent 
values of R n (although R4 drifts off by a few percent at 
K = 0.25). 

In ID the skewness amplitude is tunable and directly 
related to the amount of particle-hole symmetry break- 
ing at the microscopic level. Figs. [l6] and [l?] show how 
the microscopic particle-hole asymmetry varies as func- 
tion of K in the 2D BCSOS model. We measure several 
local quantities. Fig. [if] shows the density difference ra- 
tio, p = (ph— p v ) I '{Ph+ Pv), between local sharp hill tops, 
Ph, and local sharp valleys, pv, as seen along ID cross- 
sections of the crystal. Fig. [17] shows the density differ- 
ence phv between local sharp hill tops and local sharp 
valleys bottoms in the 2D surface and also the difference 
between sharp ridges versus sharp canyons p rc as defined 
in Fig. ^ . All three quantities vary dramatically with 
K. This demonstrates that the local skewness coupling 
constant u s k varies significantly with K. Sufficiently to 
expect a large variation in R 3 if (macroscopic) skewness 
is tunable. Since this not the case, R3 ~ —0.27 is most 
likely a universal property of the stationary state. 

The curves in Fig. |l6| and Fig. [H] have kinks at K = 0. 
These are caused by the definition of the deposition prob- 
abilities, p = min(l, exp(— AE)). The probabilities are 
temperature dependent for some configurations but con- 
stant, p — 1, for others. At K — they are being reshuf- 
fled. The curves of Fig. |l^ have similar dips at K = for 
small L, but these vanish with system size. This is an- 
other illustration of the insensitivity of the macroscopic 
length scale properties of the stationary state on the local 
properties. 



VI. CONCLUSIONS 



The most important result of this paper is that the 
amplitude ratio's R 3 = W 3 /Wl 5 and R 4 = W*/Wi 
of the stationary state moments converge to the same 
value in the KK and BCSOS model, R 3 = -0.27(1) and 
i?4 = +3.15(2). At the start of this project we expected 
that the skewness, R 3 , would be the most sensitive pa- 
rameter to test the universality of the KPZ stationary 
state properties in 2D. However, R 3 is numerically much 
more stable than the surface roughness critical exponent 



5 



a. -R3 takes the same value in the KK and BCSOS mod- 
els, and does not vary significantly in the BCSOS model 
with the temperature parameter K; in contrast to the 
strong variation in the microscopic measures of particle- 
hole asymmetry. The 2D KPZ stationary state skewness 
is universal, unlike ID where it is tunable. 

The differences in the numerical values for a in the 
KK and BCSOS model have been a puzzle for a long 
time. The corrections to FSS scaling in the KK model 
are small, compared to the MC noise and point clearly 
to the KKL value a = §. The FSS corrections to scal- 
ing in the BCSOS model are large. The stationary state 
roughness exponent converges naively but systematically 
to a smaller value a ~ 0.38. However, our FSS analysis 
shows that the data is consistent with a — | if the lead- 
ing corrections to scaling exponents are dominated by 
yi r ~ —0.6 (2) for the even moments and yi r ~ —1.7 (3) 
for the odd ones. These values are different, but in the 
same range as predicted by simple minded power count- 
ing. This is an important issue, because different values 
of a in the KK and BCSOS models, and variations with 
K in the latter, would imply non-universality of 2D KPZ 
scaling, and open up the possibility of e.g., a continuously 
varying a. Our data suggest one single KPZ fixed point 
with unique exponents. It is most likely Lassig's station- 
ary state. 

This work is supported by NSF grant DMR-9700430. 
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FIG. 1. Stationary state second moment Wi of the 2D 
Kim-Kosterlitz (KK) model as function of lattice size L 2 . The 
solid line with slope 0.8 is shown as reference. 
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FIG. 2. Finite size scaling approximants for the surface 
roughness exponent «2 of the second moment, W2 — AL a2 , 
in the stationary state of the 2D KK model. 
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FIG. 3. Finite size scaling approximants 03 (L) for the sur- 
face roughness exponent a 3 of the third moment, W3 ~ AL" 3 , 
in the stationary state of the 2D KK model. 
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FIG. 5. Finite size scaling behaviour of the skewness am- 
plitude ratio R 3 = W 3 /Wl 5 in the 2D KK model (the circles) 
and the K = BCSOS model (diamonds). 
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FIG. 4. Finite size scaling approximants on{L) for 
the surface roughness exponent Q4 of the fourth moment, 
Wi ~ AL ai , in the stationary state of the 2D KK model. 
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FIG. 6. Finite size scaling behaviour of the fourth mo- 
ment amplitude ratio R 4 = W^/Wi in the 2D KK model (the 
circles) and the K = BCSOS model (diamonds). 
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FIG. 7. Finite size scaling approximants for the surface 
roughness exponent c*2 of the second moment, W2 ^ AL" 2 , 
in the stationary state of the 2D BCSOS model. 



1.4 
1.35 

1.3 
1.25 

1.2 
1.15 

1.1 
1.05 
1 

0.95 




0.03 
1/L 



FIG. 8. Finite size scaling approximants for the surface 
roughness exponent 03 of the third moment, W3 ~ AL" 3 , in 
the stationary state of the 2D BCSOS model. 




FIG. 10. Finite size scaling approximants for the surface 
roughness exponent ai of the second moment, W2 — AL a2 , 
in the stationary state of the 2D BCSOS model as function 
of if. 
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FIG. 11. Finite size scaling approximants for the surface 
roughness exponent a.2 of the second moment in the station- 
ary state of the 2D BCSOS model at K = ±0.25 and 0. 
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FIG. 9. Finite size scaling approximants for the surface 
roughness exponent Q4 of the fourth moment, W4 m AL" 4 , 
in the stationary state of the 2D BCSOS model. 
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FIG. 12. Finite size scaling approximants for the surface 
roughness exponent Q3 of the third moment in the stationary 
state of the 2D BCSOS model at K = ±0.25 and 0. 
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FIG. 15. Finite size scaling behaviour of the fourth mo- 
ment amplitude ratio W4/W2 m the stationary state of the 
2D BCSOS model at K = ±0.25 and 0. 
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FIG. 13. Finite size scaling approximants for the surface 
roughness exponent «4 of the fourth moment in the stationary 
state of the 2D BCSOS model at K = ±0.25 and 0. 
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FIG. 14. Finite size scaling behaviour of the skewness am- 
plitude ratio i?3 = W3/W2' 5 in the stationary state of the 2D 
BCSOS model at K = ±0.25, and 0. 
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FIG. 16. The difference in the expectation value ol local 
sharp hill tops and of local sharp valley bottoms as function 
of K in the 2D BCSOS model along ID cross sections through 
the surface. The diamonds (circles) are for lattice size L = 36 
(64). 
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FIG. 17. The difference in the expectation value of sharp 
hill tops and of sharp valley bottoms, p hv , and local sharp 
ridges and canyons, p rc , as function of K in the 2D BCSOS 
model. The diamonds (circles) are for lattice size L = 36 (64). 
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FIG. 18. Definition of local sharp hill tops, valley bot- 
toms, ridges, and canyons in the 2D BCSOS model. 



10 



